Seismic exploration using non-impulsive vibratory sources activated by stationary, Gaussian codes, detecting vibrations via receivers within a wellbore and processing the results into distortion-free final records

ABSTRACT

The disclosed invention is a method of seismic exploration using non-impulsive vibratory sources activated by stationary, Gaussian codes in a wellbore environment. The method adds a distortion-free quality to the full-waveform seismic records.

CROSS-REFERENCE TO RELATED APPLICATION

This application relates to U.S. Ser. No. 316,006, filed Oct. 28, 1981, for "Seismic Exploration Using Non-Impulsive Vibratory Sources Activated By Stationary Sign-Bit Recording Gaussian Codes and Processing That Results in Distortion-Free Final Records Particularly Useful in Stratigraphic Trap Determination", assigned to the assignee of the present application now U.S. Pat. No. 4,486,866; and also relates to the following U.S. Applications filed simultaneously herewith, viz., Ser. No. 525,761 for "Seismic Exploration using Non-impulsive Vibratory Sources Activated by Stationary, Gaussian Codes, To Simulate an Impulsive, Causal Generating, Recording and Pre-processing System and Processing the Results Into Distortion-free Final Records", and Ser. No. 525,760 for "Seismic Exploration Using Non-impulsive Vibratory Sources Activated By Stationary, Gaussian Codes, and Processing the Results In Distortion-free Final Records Particularly Useful in Urban Areas".

FIELD OF THE INVENTION

This invention relates to seismic exploration in which seismic vibrations, generated by one or more non-impulsive vibratory sources, are injected into the earth, and subsequent vibrations are recorded by one or more detectors positioned within a wellbore displaced from the source. More particularly, the invention relates to the enhanced recovery of all desired seismic wave form information even though either representations of the waves injected into the earth, or representations of the subsequent vibrations of the earth, are subject to non-linear distortion. In a further aspect, the present invention has the surprising capability of providing identification of trends within the resulting records so as to indicate accumulations of oil and gas, adjacent to wellbore. In this aspect, the resulting records can be displayed as a series of amplitude-v.-time traces called "vertical seismic profiles (VSP)".

BACKGROUND OF THE INVENTION

The terms "vertical seismic profile (VSP)" used herein in its broadest sense.

In addition to the usual application of the term to describe a series of series of side-by-side traces associated with activation of single non-impulsive source, the term also is meant to describe a stacked series of traces associated with a series of non-impulsive sources which have changed lateral position with respect to the wellbore. In this regard, the processing of such data includes the steps of correlating the input and received date; separating the up- and down-going waves after correlation and then displaying the up- and down-going waves as a function of time for the purpose of structural identification to allow orderly development of the oil and gas accumulations adjacent the wellbore. As to examples of processing related to the two last-mentioned steps, see "COMPUTER PROCESSING OF VERTICAL SEISMIC PROFILE DATA", M. W. Lee and A. H. Batch, GEOPHYSICS, V. 48, No. 3, 1983.

Since shape and persistence of the waveforms as a function of time are most important characteristics of accurate vertical seismic profiles, where distortion resides within the vibratory system, the resulting traces have been found to be most difficult to interpret. That is to say, accurate VSP interpretation requires that the compared traces be essentially distortion-free, after processing has occurred.

The term "distortion-free, after processing" to describe the final records relates to the fact that there should be no undesired or unknown, linear or non-linear change in the down-going waveforms either of the recorded field signals or of the final data after processing has been completed.

Prior to my invention, such distortion could be brought about (i) by the use of techniques in the field as described above as well as hereinafter wherein harmonics of the fundamental dissipate its energy, or (ii) accidentally, where for example, vibrator pad and the earth do not linearly interface.

Doty et al., U.S. Pat. No. 2,688,124, issued Aug. 31, 1954, for "Method and Apparatus for Determining Travel Time of Signals" describes the well-known VIBROSEIS® system of Continental Oil Company. In such a system, seismic waves are generated by mechanical vibrators on the earth's surface. Each of the vibrators is firmly anchored to the earth by the combined weight of the source. Peak forces in the neighborhood of 10 to 20 tons (and up to 36 tons) can be developed by the rapid, non-explosive interaction of the base-plate and piston system of each vibrator. Consequently, the weight of each vibrator is proportionally large to maintain the desired, continuous vibrator-earth contact during operations. The waves sent into the earth consist of long, sinusoidal wave, trains of predetermined frequency and time duration characteristics rather than the much sharper wave impulses sent into the earth by the explosive sources used prior to the VIBROSEIS® system, or by "weight drop" methods including those provided by various impulse-coded systems, e.g., the so-called "pulse-coded" techniques.

There may be some confusion as to the differences of the signals produced by the VIBROSEIS® system and those produced by impulsive sources such as provided by exploding dynamite, exploding mixtures of propane and air, or by "weight drop" methods, including pulse-coded techniques.

It is well known that the capacity of any signal (including seismic signals) to carry information can be measured in a manner analogous to determining the volume of a container. Since volume is the product of height times width times length; similarly, information capacity of a signal is related to a product of amplitude, frequency bandwidth and the length of the signal.

Dynamite as a seismic energy source produces an input signal having considerable amplitude (height) and bandwidth, but has very short length. On the other hand, "non-impulsive" vibrations generated in the manner of a VIBROSEIS® system have limited amplitude, but such is compensated for by the long length of the input signal and a faithful, continuous reproduction of the control signal over the frequency spectrum of interest. That is to say, in the VIBROSEIS® system, the amplitude and phase spectra are carefully and continuously controlled so that the resulting energy spectra changes smoothly as a function of time. Thus, a smoothly varying output of desired frequency and duration characteristics is provided in contradistinction to the binary-coded (ON-OFF) square wave output generated by pulse-coded methods in which the energy-per-blow is substantially constant and cannot be so controlled.

A further essential part of the VIBROSEIS® system lies in the processing of the received data to produce records that tend to show short pulses representing reflections from subsurface interfaces. Such responses are provided by cross-correlating the recorded representation of the vibratory waves sent into the ground with the recorded representation of the waves received subsequently.

The use of cross-correlations, as taught by Doty et al. and many others since, has now become so well known in vibratory seismology that it will be presumed to be well known in the following parts of the present specification; and the description will concern itself only with differences from the prior art.

Erich, U.S. Pat. No. 4,234,053, for "Seismic Exploration Method Using a Rotating Eccentric Weight Seismic Source", describes an exploration method in which a rotating eccentric weight source is used (as a power impactor) to transmit a coded, non-Gaussian impulse input signal into the earth on a substantially constant energy-per-blow basis. An improved representation of the pulsed input signal is correlated with the raw seismic data to provide the field record of interest. But since the impulsive source is also only discontinuously coupled to the earth (i) the interaction of the mass of the eccentric weight source with the spring constant of the earth produces an output dominated by low frequency components and (ii) the pulse shape of the output can vary non-linearly with time. Hence, such a system is limited to those uses where a conventional VIBROSEIS® system cannot be employed.

Multi-array use of such sources is likewise limited.

Another relevant patent is that of Crook et al., U.S. Pat. No. 3,264,606, issued Aug. 2, 1966, for "Method and Apparatus for Continuous Wave Seismic Prospecting", which teaches driving of vibratory sources (in conjunction with conventional full-wave recording equipment) with pseudo-random codes which, although differing in detail from the preferred codes prescribed here, does share the desirable generic property of "a code sequence which may be represented as a reference time series having a unique auto-correlation function comprising a single major lobe having no side lobes of greater amplitude than the side lobes of the auto-correlation function of statistically unrelated noise components of the composite signal detected at said detecting location" (column 13, lines 32-44).

In my U.S. Pat. No. 4,346,461 for "Seismic Exploration Using Vibratory Sources, Sign-Bit Recording, and Processing That Maximizes the Obtained Subsurface Information", issued Aug. 24, 1982 and assigned to the assignee of the present application, I describe a non-impulsive vibratory system that uses a class of vibrator signals best characterized as Gaussian, zero-mean, and stationary, in conjunction with sign recording of both the injected and received vibrations at the sources and receivers. The stated advantages relate to the channel-capacity economy of sign-bit recording (at both the sources and receivers), and to the distortion-free quality of the final processed records.

I have now discovered that use of the above class of vibrator signals not only does not sacrifice information in the final processed records even though the data is collected by sign-bit recording methods, but such type of vibrator signals also favorably impacts seismic recording and processing operations within wellbore environments where full-wave recording is contemplated.

Aside from the above, a paper of A. B. Cunningham, Geophysics, December 1979, Vol. 44, No. 12, pages 1901 et seq., for "Some Alternate Vibrator Signals", works out in mathematical detail expected types of cross-correlation functions from various types of vibrator sweeps, including certain types of pseudo-random sweeps, but not in the context used herein.

SUMMARY OF THE INVENTION

The present invention uses a class of vibrator signals best characterized as Gaussian, zero-mean, and stationary, in conjunction with conventional recording techniques in a wellbore environment. In this aspect of the present invention, the distortion-free nature of the final displays results even though the vibratory system may not be distortion-free during the injection of the above-identified vibrations into the earth adjacent to the wellbore and recording of a full-wave representation of the injected signal. The vibrations are associated with a zero-mean, stationary, Gaussian code used to drive the vibrator. Received signals are next recorded at one or more receivers positioned at known locations within the wellbore. Full-wave preservation methods of recording are used. Thereafter, the recorded signals are processed using conventional cross-correlation techniques including displays of the waves on a side-by-side basis normalized to sourcepoint location, for example. Absolute scale of each received trace is maintained.

Since the final records are surprisingly distortion-free, they have a superior capability of, say, indicating oil and gas accumulations associated with the surveyed earth formation adjacent the wellbore.

These and other features and advantages of the present invention will become evident to exploration seismologists in light of the following specification, including the appended drawings briefly described here.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a side elevation of an earth formation penetrated by a wellbore and shows major pieces of apparatus used in a survey system for carrying out the present invention, the survey system including a vibratory source adapted to be driven by a Gaussian, zero-mean, stationary code, a series of receivers located in the wellbore to receive the propagated portions of the injected vibrations, and processing and recording means for the generated and recorded signals.

FIG. 2 shows a sequence of auto-correlations, first that of one random sweep, then a sum of the auto-correlations of two different random sweeps, then a sum of the auto-correlations of 4, then 8, then 16 different random sweeps.

FIG. 3 is a schematic diagram in block form of the system of FIG. 1 showing paths along which signal information flows, in accordance with the present invention.

DETAILED DESCRIPTION OF THE INVENTION

Refer now to FIG. 1, 10 illustrates the survey system of the present invention in a wellbore environment. System 10 includes vibratory source 11 positioned adjacent to a wellbore 12, say at a known sourcepoint location (SP). Within the wellbore 12 are one or more detectors 13. In operation, vibratory source 11 injects seismic vibrations in a continuous and non-impulsive manner into the earth formation 14. Amplitude-versus-time characteristics of the vibrations of the source 11 vary over a predetermined time duration and frequency spectrum, under control of a control signal generated by sweep generator 15. The sweep control signal fed to vibratory source 11 is also fed through a separate information channel to correlating and recording unit 17.

In FIG. 1 the vibratory wave going into earth formation 14 are represented by ray paths 20, the waves traveling obliquely downward toward a pair of stratal interfaces 21, 22 where they undergo partial reflection and travel back upwards along paths 23, again obliquely, toward the detector 13. In addition to the detecting waves traveling upwardly from interfaces 21, 22 (called "up-going waves"), the detector 13 within wellbore 12 also collects waves that travel downwardly ("down-going waves"), as along path 24. While the down-going waves traveling along path 24 are useful in determining the velocity characteristics of the formation situation between the sourcepoint (SP) and the detector 13, such waves can interfere with upcoming waves arriving at the detector 13 and must be identified and subsequently separated therefrom. As discussed below, such identification and separation can conveniently occur after cross-correlation of the full-wave received and input signals within processing and recording unit 17, has occurred.

While only one symbolic detector is shown in the wellbore 12 of FIG. 1, in actual current field practice, several detectors may be employed. From detector 13, the final operations represented in FIG. 1 are transmissions of full-waveform data to the correlating and recording unit 17. In the unit 17, cross-correlations are made of the full-wave received data and the injected vibratory sweeps.

In this regard, when the resulting waveforms were cross-correlated, cross-correlations from injected vibrations associated with Gaussian, zero-mean, stationary driving codes are superior to those associated with the chirp signal. But the mere appearance of superiority would not convince seismologists in general that the desirable information is due to the above-mentioned drive code especially where the vibratory system, itself, is not distortion-free.

But it is now possible to state that for currently practiced seismic exploration in a wellbore environment, if the specified types of driving functions are used, as taught in the present invention, the resulting seismic information will be surpisingly distortion-free.

The physical steps of the present invention, when properly combined, make possible the utilization of the theoretically provable, maximum amount of information that can be recovered from the correlated waveforms.

That it should be possilbe to better recover distortion-free information from correlations of a stationary, zero-mean, Gaussian function was obliquely implied in theoretical results obtained by Bussgang. These results are concerned with noise for the purpose of data communication and, to my knowledge, have remained relatively unknown among exploration geophysicists. Their geophysical implications require an unobvious twist of thought in that it is necessary to think of Bussgang's "noise" as a possible seismic "signal", just the opposite from what would then be called the "noise" in the seismic context. [Chapters 5 and 6, e.g., of An Introduction to Statistical Communication Theory, J. B. Thomas, John Wiley & Sons (1969), are seen to express "noise" and "signal" in such conventional terms.]

The results of Bussgang that bear upon the present invention are equations given here in notation different from that originally used notation, but more like that used in current theoretical seismology.

Bussgang's result can be stated as follows:

Let

W, Y be stationary, jointly Gaussian, zero-mean, random processes (continuous in nature or sampled representations of such continuous processes);

X(.,.) be the cross-correlation operator; and

k(.) be a characteristic of any non-linear function (i.e., k is a fixed mapping function or table that converts input values to output values but not necessarily one-to-one).

Then

    X(W,Y)=a; X[W,k(Y)]=b; X[k(W), Y]                          (2)

where a and b are constants depending only on k and the sequence affected by k.

In words, equation (2) states that given two input signals (W and Y) with certain properties, their cross-correlation X(W,Y) is unaffected, except to a scale factor, by any distorting mechanics that either one might suffer.

In such case, X(W,Y) represents the correlated output of a true amplitude recording system (i.e., without non-linear distortion); a X[W,k(Y)] represents the correlated output of a system in which the data undergoes non-linear distortion as it comes through the geophone, and then is correlated with the sweep signal; and b X[k(W,Y] represents the correlated output of a system in which the sweep undergoes non-linear distortion, and then is correlated with the unclipped output of the geophone.

Interpreted in this fashion, we can distinguish the two cases where:

(i) Y, the recorded signal, is subjected to the non-linearity, k, before correlation; and

(ii) W, the source signal, is subjected to k before correlation.

These results apply directly to seismic prospecting where W is the motion of a vibratory source and Y is the consequent motion of the earth picked up by a geophone of a series of geophones.

In this case it is well known that W and Y are related through the usual convolutional equation:

    Y=E*W+N                                                    (3)

where E is the impulse response of the earth and N represents any uncorrelated noise present at the geophone.

Interpreted in this fashion, the extension of Bussgang's result to seismic prospecting implies that if the sweep and recorded data have the right properties, non-linear operations produce no distortion; but in case (i) above there is a loss of absolute trace amplitude information. Absolute amplitude data is lost because Bussgang's extension applies to normalized cross-correlations only. The loss is the same as that which occurs when data are trace equalized using a single gate containing the whole trace. It is not equivalent to an AGC operation. In case (ii), the correlated output is undistorted and, since the scaling constant is only a function of the particular non-linearity and the particular sweep statistics, if these are held constant from profile-to-profile, then true scaling both for the traces within a profile and for the traces from different profiles, is assured.

OPERATIONAL STEPS OF THE INVENTION

With the above as a guide, briefly, the operational steps of method of the present invention may now be specified using vibratory sweeps and conventional full-wave recordings, so as to produce final cross-correlation traces that are surprisingly distortion-free and so are superior to conventional cross-correlation traces. In accordance with method aspects, a preferred embodiment of the present invention includes the steps of:

A. Injecting seismic vibrations into the earth by driving at least one vibratory source using a control code of preselected Gaussian, zero-mean, stationary characteristics;

B. Recording the full-wave representations of the injected seismic vibrations and received full-wave vibrations that have propagated through the earth from the vibratory source(s) to the receiver(s) within an adjacent wellbore; and

C. Cross-correlating the full-wave representations of the injected vibrations with the received signals.

The result of such cross-correlation, without more, is the seismic record of interest that is surprisingly useful in oil and gas exploration.

The above three steps, performed in combination, give the result never achieved before: the final cross-correlation traces that are superior to cross-correlation traces that could have been made heretofore using an impulsive or a chirp source signal.

Briefly, superiority of the method of the present invention results from the fact that the use of the Gaussian, stationary vibrations invert the signal collection and processing requirements for producing distortion-free final correlated records. Thus, where the prior art performs any form of non-linear operation on the source or receiver signals, distortion would occur. But note that in my method, all such non-linear operations are performed in an inverse domain (the "noise" domain) and hence do not distort the essential nature of the final seismic records. Hence, the present invention can operate with superiority in near surface operations involving vertical seismic profiling operations.

Within the above-identified area of usage, the step in the present invention, related to the injection of vibrations into the earth using a Gaussian, zero-mean, stationary driving code for the vibrator, plays a very important role.

Hence, it is appropriate therefore to go into some detail concerning the nature of such vibrations and especially the control codes used in conjunction with sweep generator 15 to drive the vibratory source 11 of the system shown in FIG. 1.

In the art of driving vibratory sources, extensive use has already been made of activating codes embodied on magnetic tape, so it is unnecessary to discuss in detail how the activating codes are used to drive the vibratory source 11. It should be commented, however, that sufficiently detailed equipment is now presently available for this purpose. For example, in Payton, Waters and Goupilland's U.S. Pat. No. 4,168,485 for "Simultaneous Use of Pseudo-random Control Signals in Vibrational Exploration Methods", a digital sweep generator is described that uses a digital control code to produce the desired control signal. Essentially, the digital control signal is placed in the memory of a digital computer as digital words which can be shifted out, in seriation, to a digital-to-analog converter. The converter, in turn, generates an alternating signal to drive the vibratory source, such a vibratory source 11. The output signal from the D/A converter is initially in the form of a series of stepped electrical functions in which amplitude changes, in steps, as a function on incremental time change. But after smoothing via a filter, the stepped signal becomes a more conventionally looking alternating control signal those amplitude varies smoothly with time, to provide the desired output signal.

It should also be noted that the art has further progressed where it is no longer necessary to use main-frame type computers on-board the field equipment truck. Advances in both areas of software and hardware technology now allow the explorationist a choice of numerous mini- and micro-computers for this purpose. Additionally, he is no longer required to generate the activating codes at a remote computer center and transport them on tape to the field location. Codes can now be generated in real time at the field site using the above-identified equipment and it seems possible that in the future this will become the preferred way to generate them.

The discussion here will concern itself with some of the possible options in preparing a code particularly suited for the practice of this invention; a code that results in the injection of an associated signal into the earth. Its duration can be longer than that of conventional "chirp" signals. For discussion purposes, 32 seconds is a possible time interval. The magnetic tape devices used in geophysical prospecting in urban areas have various sampling intervals. Two milliseconds is a possible choice. Thirty-two seconds of code at a two-millisecond sampling rate require 16,000 random numbers. These may be generated in several ways. One way is to use, in succession, two of the well-known subroutines in the IBM Scientific Subroutine Package: RANDU, to generate a set of uniformly-random numbers, then GAUSS, to convert those numbers into a set with Gaussian distribution. Other ways of generating suitable sets of numbers have been discussed by Lewis in the IBM Systems Journal No. 2 (1969) and by Knuth in his book "Seminumerical Algorithms" (Vol. 2 of the "Art of Computer Programming", Addison Wesley Pub. Co.). However it is done, the desired objective for purposes of this invention is a set of random numbers with a zero-mean, Gaussian, stationary vibratory signal.

Another way is to employ a controlled code of preselected characteristics, but which is modified in a fast, feedback mode as a function of the output vibrational signals. Result: irrespective of the characteristics of the control code, the total vibratory system interacts so as to provide the desired output signal.

A set of random numbers, 500 per second, may contain, in frequency terms, components up to 250 Hertz, and such frequencies are higher than those usually considered useful in current field vibrator work. So a possible next step is to put the set of random numbers through a digital band-pass filter to restrict the frequency content to a band, say, between 10 and 80 Hz. Many phase characteristics are possible.

The last condition that the code should satisfy in order that the entire method of the present invention will correspond as nearly as practicable to the ideal conditions indicated by Bussgang's results, is the condition of "stationarity". Statistical stationarity, by definition, is the maintenance of the same statistical characteristics throughout the time interval of interest. A signal generated as described above has inherent stationarity, e.g. as long as program control parameters are not varied during the time it is activating the vibrator(s). Stationarity of both the source and receiver data is approached as closely as practicable by (1) beginning the actuation of the vibrator(s) long enough before the starting of the recording for transients associated with the impulse response of the earth to die down substantially (e.g., 6 seconds), and (2) continuing the non-impulsive actuation of the vibrators, such as 11, throughout the long-time interval conventional for the vibrations to propagate through the earth to the receiver, such as detector 13 within wellbore 12.

It was stated above, without the benefit of any supporting discussion, that the vibratory source(s) should be "substantially linearly-responding" (if controlled by a Gaussian, zero-mean, stationary driving code, for example), or the entire vibratory system can interact within itself, to provide the desired vibratory output. These conditions for properly carrying out the present invention may require changes from some of the present field practice of vibratory sources.

EXAMPLE

A field test is illustrative of the method of the present invention. In the test, the receiver nearest the shot is offset from the source. Spacing between receiver is 30 feet.

Structure undergoing survey in the present example is conventional structural situation.

In one instance, comparative data is generated using a conventional "chirp" signal of about 16 seconds duration of the form: ##EQU1## where the signal consists of a sinusoidal function of gradually and uniformly increasing frequency, starting at w_(o) =2πa(f_(o) =a) and ending at w_(f) =2π(a+bT). Thus, a sweep going from 10 to 100 Hz in 16 seconds has values of a=10 and b=5.6. The purpose of the function A(t) is to have a multiplying factor in the sinusoidal function so that the sweep does not start or stop abruptly in amplitude.

The spectrum of the system input is assumed to be completely white so the shape of the amplitude spectrum of the received signal is the shape of the earth signal.

In the other instance, the driving code for the vibrator is Gaussian, zero-mean, and stationary. Duration: 16 seconds long. Sample rate: 2 milliseconds. The sweep is generated by passing a random number sequence generated with RANDU and GAUSS through a zero-phase filter. The amplitude response of the filter increases linearly from 0 Hz to 125 Hz at a rate of 12 db/octave. Above 125 Hz, the filter amplitude response decreases linearly at a rate of 72 db/octave.

Both sets of data are correlatable. I.e., to say a full amplitude representation of the conventional "chirp" signal is cross-correlatable with the full-wave signal received at each receiver group. Similarly, the full-wave representation of the Gaussian, zero-mean, stationary sweep signal is cross-correlatable with the full-wave receiver signal of each receiver.

It is evident from the data developed that the method of the present invention is far superior to processes using the conventional "chirp" method. Better resolution of adjacent waveform changes (in the traces provided by the method of the present invention) is the main advantage. Hence there is strong evidence that the present method provides final records that are superior to final processed records obtained using a conventional "chirp" signal, especially for indicating the presence of structure in a wellbore environment.

From this point on, the discussion will be concerned with refinements and variations not necessary to make the basic method operable, but may be desirable for optimum performance.

Sometimes in vibrator seismology, in order to increase the effective signal-to-noise ratio, it is desirable to make repetitive sweeps and add the records from the successive repetitions expecting that desired signals will be reinforced and the noise will tend to cancel itself out. When Gaussian, zero-mean, stationary sweeps are being used, and such repetitions are made, it is desirable to use a different time sequence of vibrations for each repetition. FIG. 2 illustrates this, using auto-correlated random sweeps recorded as traces 36a, 36b, 36c, 36d, and 36e. Trace 36a is the auto-correlation from a single sweep. The other traces 36b, 36c . . . 36e are the sums of the auto-correlations of successively 2, 4, 8, and 16 different sweeps. (The vertical lines are 100 millisecond timing lines.) The events seen in addition to the main peak are due to correlation noise. Since the correlation noise of any one sweep is different from that of any other sweep, the level of this noise decreases as more sweeps are summed.

Some of the variations and alternatives in the method of the present invention are illustrated in FIG. 3. Alternate data paths are shown along which signal information may flow during practice of the method. The original Gaussian code from the sweep generator 15 may flow through filter 37 on its way to the vibrator 11, or it may not, depending on the state of switch 38a. The same code may or may not flow through filter 39, depending on switch 38b. Subsequently, the signal passes to the correlator and recorder unit 18. Representations of the full-wave received signals are seen to pass into the recorder unit 18. Within unit 18, the multi-bit correlation functions go directly to form the final record.

The invention is not limited to the specific embodiments described above, because variations will be readily apparent to those skilled in the art. Thus, the invention is to be given the broadest possible interpretation within the terms of the following claims. 

What is claimed is:
 1. A method of seismic surveying using non-impulsive vibratory sources that recovers all the essential information of full-waveform recording in a distortion-free manner in a wellbore penetrating an earth formation, which comprises:A. injecting non-impulsive, seismic vibrations into the earth by activating at least one vibratory source with a driving code of preselected Gaussian, zero-mean, stationary characteristics; B. recording the seismic vibrations that have propagated through the earth from said at least one vibratory source to at least one receiver within said wellbore; and C. cross-correlating a full-waveform preserved representation of said injected vibrations with received seismic vibrations; the result of said cross-correlation being a seismic survey record that is essentially distortion-free.
 2. The method of claim 1, in which the time duration of said seismic vibrations is continuous throughout the time interval in which said propagated seismic vibrations are received, and essentially is equal to that for a conventional "chirp" signal.
 3. The method of claim 1, in which when Steps A to C are repeated, a different time sequence of Gaussian, zero-mean, stationary driving codes is used for each time said steps are repeated.
 4. The method of claim 1, in which said driving code of Gaussian, zero-mean, stationary character is filtered before said code is used to drive said vibratory source, using a filter selected from the class consisting of minimum-phase band-pass filters and zero-phase bandpass filters.
 5. The method of claim 1, in which said driving code of Gaussian, zero-mean, stationary character is filtered before said code is used to drive said vibratory source and also before said code is recorded, through a filter selected from the class consisting of minimum-phase band-pass filters and zero-phase band-pass filters.
 6. A method of seismic surveying in an earth formation penetrated by a wellbore, comprising:A. generating a Gaussian, zero-mean, and stationary code; B. filtering said code through a filter selected from the class consisting of minimum-phase band-pass and zero-phase band-pass; C. driving at least one substantially linearly-responding seismic vibratory source with the resultant filtered code to inject vibratory signals into the earth for a time interval longer than the recording interval of the resultant seismic waves; D. recording the full-wave seismic signals detected by at least one receiver positioned in said wellbore to detect said resultant seismic waves; E. recording a waveform-preserved representation of said Gaussian, zero-mean, stationary code without filtering thereof; and F. cross-correlating said waveform-preserved representation of said code with the full-wave detected seismic signals; the result of said cross-correlating being a seismic record that is essentially distortion-free.
 7. The method of claim 6, in which Steps A to F are repeated, and a different time sequence of Gaussian, zero-mean, stationary seismic driving codes is used for each time said steps are repeated.
 8. A method of seismic surveying using at least one, non-impulsive vibratory source that recovers all essential information of full-waveform recording, which comprises:A. injecting non-impulsive, seismic vibrations into the earth by activating at least one vibratory source with a driving code of Gaussian, zero-mean, stationary characteristics; B. recording a full-wave representation of said injected vibrations; C. recording a waveform-preserved representation of the seismic vibrations that have propagated through the earth from said at least one vibratory source to at least one receiver in said wellbore in an up-going direction; and D. cross-correlating said full-wave representation of said injected vibrations with said full-waveform-preserved representation of said propagated seismic vibrations recorded by said receiver; the result of said cross-correlation being a seismic survey record essentially distortion-free and particularly useful in structural determination in said surveyed area adjacent said wellbore. 